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Abstract. We study the correlation and response dynamics of trap models of glassy 
dynamics, considering observables that only partially decorrelate with every jump. 
This is inspired by recent work on a microscopic realization of such models, which found 
strikingly simple linear out-of-equilibrium fluctuation-dissipation relations in the limit 
of slow decorrelation. For the Barrat and Mezard model with its entropic barriers we 
obtain exact results at zero temperature T for arbitrary decorrelation factor k. These 
are then extended to nonzero T, where the qualitative scaling behaviour and all scaling 
exponents can still be found analytically. Unexpectedly, the choice of transition rates 
(Glauber versus Metropolis) affects not just prefactors but also some exponents. In the 
limit of slow decorrelation even complete scaling functions are accessible in closed form. 
The results show that slowly decorrelating observables detect persistently slow out-of- 
equilibrium dynamics, as opposed to intermittent behaviour punctuated by excursions 
into fast, effectively equilibrated states. 



1. Introduction 

Trap models [1112113111130111013 UHl Ulj have been recognized in recent years as 
powerful models of glassy, non-equilibrium dynamics. They describe the motion of a 
system through its phase space, simplified to a picture of thermally activated hopping 
in a landscape of traps of energy E. The simplest case is that of mean-field trap models, 
where all traps are taken as mutually accessible and the rate for a transition between 
two traps depends only on their energies. An alternative motivation for considering 
trap models is provided by the problem of diffusion in disordered media; here the traps 
are located in a physical space of low dimension (say, d = 1,2 or 3) and the spatial 
organization of the traps has to be accounted for [121 CSl IHl HHj- 

Our motivation in this paper arises from a recent interesting study of a microscopic 
realization of trap models on the basis of the number partitioning problem P^ ITT]. Each 
partition of N numbers between two piles can be associated with the state of a system of 
Ising spins Sj = ±1, and an energy function is then defined to measure how far from 
optimal the partition is. The dynamics considered by Junier and Bertin in jT2j is that at 
every step K spins are selected and have their values randomized; the new state is then 
accepted with a probability given by a Metropolis acceptance factor. Junier and Bertin 
argue that this system evolves in a manner analogous to the trap model considered by 
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Barrat and Mezard [2]: a new state or trap is selected essentially at random at every step, 
and accepted according to the Metropolis (or Glauber, see below) probability. To justify 
the assumption that each new state is effectively random one requires ^ 1; for finite 
K there is a crossover time, exponentially large in K, beyond which this simplification no 
longer applies [TJj. In contrast to the usual assumption made in trap models, however, 
an observable such as the (randomly staggered) magnetization does not decorrelate fully 
with every transition between states. Instead, it decorrelates by a factor k = 1 — K/N, 
since only K out of the spins are updated. Junier and Bertin showed that the 
consequences of this are rather profound, and most interesting in the limit of slow 
decorrelation n —>■ 1 {K <^ A^). In particular, they found that fluctuation-dissipation 
(FD) relations between response and correlation functions had the simple straight-line 
form expected for systems with well-defined (effective) temperatures [TSt ITOt l2n| l2T| . 
The corresponding temperature was equal to the bath temperature T down to half the 
glass transition temperature T^; for lower T, it remained pinned to T^s = Tg/2. This 
unusual transition between apparent equilibrium and non-equilibrium dynamics within 
the glass phase was interpreted as due to a change from activated to entropic slowing 
down. 

Our aim in this paper is to complement the work of [^, which used mainly 
numerical simulation and simple scaling estimates, with an analytical study of the 
dynamics of slowly decorrelating observables in trap models. This will allow us, for 
example, to verify and refine estimates of scaling exponents in [T7j, but also yield more 
detailed insights into the nature of the dynamics. In our approach it is also a simple 
matter to compare different choices of the transition rates (Metropolis and Glauber) 
and we will see that this has some unexpected consequences, with effects not just on 
prefactors but also scaling exponents. 

To allow a direct comparison with the work of we focus mainly on the Barrat and 
Mezard model itself; the extension of this to slowly decorrelating observables is defined 
in Sec. |21 General expressions for correlation and response functions are then derived 
in Sec. 01 and simplified in Sec. EJfor the interesting scaling regime of long times. As a 
prelude to the analysis of slowly decorrelating observables proper, we consider in Sec. El 
first the standard, fully decorrelating case k = at nonzero temperature, extending 
previous analytical results for the FD behaviour in the limit T — ^I]. Sec. El then 
looks at arbitrary decorrelation factors k, but first at zero temperature, where closed 
form results can still be obtained. The most general case of k > and T > is then 
analysed in Sec. [71 where we also provide results from numerical solutions of the integral 
equations for the scaling functions. In Sec. |Hl finally, we extend the analysis to more 
general trap models with a multiplicative dependence of the transition rates on the 
applied field jTUI; this covers in particular Bouchaud's original trap model P with its 
purely activated dynamics. Sec. IHl summarizes the results and discusses the dynamics 
of slowly decorrelating observables in a wider context. 



Trap models with slowly decorrelating observables 



3 



2. Model definition 

Motivated by the arguments discussed in the introduction, we consider a modified Barrat 
and Mezard model j2] where the magnetizations before and after a jump between two 
traps are at least partially correlated. Transitions can still take place between arbitrary 
trap energies, but are more likely between traps with magnetizations that are sufficiently 
close to each other. 

The dynamics of the system is described by the probability P{E, m, t) that the 
system is in a trap with energy E and magnetization m at time t. Of course m is 
in principle a generic observable, but we continue to use the term magnetization for 
definiteness. We assume that, starting from a trap with energy E' and magnetization m', 
a jump to a new trap with energy E and magnetization m is attempted with probability 
p{E) p{m\m') . The first factor is simply the density of states of trap energies, refiecting 
the assumption that transitions between arbitrary energy levels are possible. The factor 
p(m|m'), on the other hand, can incorporate various degrees of correlation between m' 
and m; if m is typically close to m', then decorrelation of m over time is slow. The 
factorization p[E) p{m\m') assumes that E and m are selected independently, which in 
a spin system is a reasonable assumption if m is some staggered magnetization that is 
uncorrelated with the energy. In a continuous-time description, the time evolution of 
P{E,m,t) is then governed by the master equation 
d 

—P{E,m,t) = -T{E,Tn)P{E,Tn,t) + 

+ JdE'dm' p{E)p{m\m')w{E - E' -hm + hm')P{E',m',t) (1) 

Here w{AE) is the probability with which a proposed transition is accepted. In 
writing the energy change AE = E — hm — {E' — hm') on which this depends we 
have allowed for the presence of a field h thermodynamically conjugate to m, in order 
to be able to deduce response properties of m. We consider primarily the Glauber 
form w{AE) = 1/[1 + exp{pAE)] of the acceptance probability that was used in 
previous theoretical studies [3 HH 123, but also compare with the Metropolis choice, 
w{AE) = 1 for AE < and w{AE) = exp{—l3AE) for AE > 0, which is common in 
simulations Here and throughout (3 = 1/T denotes the inverse temperature. In 

the first term on the right of (JJ, finally, we have defined 

r(E,m) = JdE' dm' p{E')p{m'\m)w{E' - E - hm' + hm) (2) 

which can be thought of as the total rate of leaving a trap with energy E and 
magnetization m. 

Our dynamics should be capable of describing thermodynamic equilibrium at 
high temperatures, and thus obey detailed balance. In line with our assumption of 
independence of energies and magnetizations in an attempted jump, let us assume the 
overall density of states also factorizes into p{E)p{m). It is then easy to check that 
the dynamics (P) obeys detailed balance as long as p{m'\m)p{m) is symmetric under 
interchange of m' and m. To be specific, we take p{m) to be a Gaussian with zero 
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mean and unit variance, and p{m'\m) a Gaussian with mean urn and variance 1 — k^. 
(Since we will only need the second order statistics of m and m', this assumption in 
fact constitutes no loss of generality.) The original Barrat and Mezard model then 
corresponds to Glauber dynamics with k = 0, while in the opposite limit k = 1 the 
magnetization remains frozen to its initial value. For intermediate values, m decorrelates 
by a factor k with each jump, and we expect the interesting behaviour seen by Junier 
and Bertin [T7j to occur in the limit of slow decorrelation k, 1 (but keeping n < 1). 



3. General expressions for correlation and response 

For the calculation of the correlation function of m, the field can be set to zero. The 
master equation then simplifies to 

^P{E, m, t) = -T{E)P{E, m, t) + J dE' dm' p{E)p{m\m')w{E - E')P{E', m', t) (3) 
with 

r(^) = J dE' p{E')w{E' - E) (4) 

The propagator P{E, 171,1 — t„|£'„,mw) obeys this equation with the initial condition 
P{E, m, 0\E^, TTiw) = S{E — E^)6{m — m„) and gives the correlation function as 

C{t,t^) = JdE dmdE^dm^mm^ P{E,m,t — t^\E^,m„)p{m^)P{E^,t^) (5) 

Here we have used that, for h = 0, P{E„,mvj,t„) = p{m„)P{E^,t„); this is true as 
long as the initial condition has the same structure, e.g. for a quench from thermal 
equilibrium at high temperature. From Q, the function p{E,t — t„\E^,m^) = 
^dmmP{E,m,t — t^\E^,m^) obeys, using ^dmmp{m\m') = Km', 

^p{E,t\E„,m^) = -T{E)p{E,t\E„,m^) + K J dE' p{E)w{E - E')p{E' ,t\E^,m^) (6) 

Because this equation is linear, the factor m„ from the initial condition 
p{E,0\E^,m„) = m„S{E — E„) pulls through and we can write p{E,t\E^,m^) = 
m^p{E,t\E^,l) = m^p{E,t\E^), dropping the constant argument 1. The reduced 
magnetization p{E,t\E„) then obeys 

—p{E,t\E^) = ~T{E)p{E,t\E^) + K JdE' p{E)wiE ~ E')p{E',t\E^) (7) 

with p{E,0\E^) = 6{E — E„). For our purposes more useful, however, is the 
corresponding backward equation, 

—fi{E,t\E^) = -T{E^)p{E,t\E„) + K JdE' p{E,t\E')p{E')w{E' - E^) (8) 
This is because the correlation function can be written as 



with 



C{t, U) = J dE^ M{t - U\E^)P{E^, U) (9) 
M{t-U\E^) = jdEp{E,t-U\E^) (10) 
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By integrating (jSI) we see that this obeys 

-M(t|E^) = -TiE„)Mit\E„) + KjdE' M{t\E')p{E')w{E' - E^) (11) 

with M(0|i?w) = 1; the forward equation ((7j), on the other hand, would yield 
an expression for dM/dt which still involves p{E,t\E^). The physical meaning of 
M{t — t^\E^) is as follows: if we start in a state with energy E^ and magnetization m^, 
then m^M(t — t^\E^) is the average magnetization a time t — t„ later. 
For the case k = the solution of (jHl) and (fTT|) is trivial, 

piE, t\E^) = 5{E - E^) exp(-r(E^)t) (12) 
M{t\E^) = exp(-r(E„)t) (13) 

and one retrieves the standard result for the hopping correlation function ITI H ITT ] I22j . 
For K > 0, the Laplace transform (LT) of (jH)) can be easier to work with; with s conjugate 
to t, this reads 

[s + T{E)]p{E, s\E^) - 5{E - E^) = k j dE' p{E)w{E - E')p{E', s\E„) (14) 

By the same reasoning that lead to ©, the response to a field impulse of amphtude 
h and duration At, applied at time t^, can be written as 

hAt R{t, t„) = JdEdmm M{t - t„\E)P{E, m, t„ + At) (15) 

The change in P{E,m,tyj + At) from its value without the field (the latter does not 
contribute to R because it is symmetric in m) is 

AP{E,m,t„ + At) = At JdE^dm„ X (16) 

[p{E)p{m\Tn„)Aw{E - E„ - hm + hm„)p{m„)P{E„, t„) 

—p{E^)p{m^\m)Aw{E^ — E — hm„ + hm)p{m)P{E, t„)] (17) 

where Aw{E — E^ — hm + hm„) = w{E — E^ — hm + hm^) — w{E — E^) is the change 
of the acceptance probability caused by the field. Expanding this to linear order in h as 
Aw{E — E„ — hm + hm„) = h{m^ — m)w'{E — E^) and carrying out the integrations 
over m and gives 

R{t,t^) = (1- k) jdEdE^M{t-t^\E)x 

X [-w'{E - E^)p{E)P{E^, t„) - w'{E^ - E)p{E^)P{E, t^)] (18) 

Note that w'{-) is negative since the acceptance probability w{-) decreases with 
increasing energy change; thus both terms in the expression for the response are positive. 
Explicitly, we have for the Glauber case —w'{AE) = P exp{PAE) /[I + exp{PAE)]'^ and 
for Metropolis -w'{AE) = Q{AE)P exp{-pAE), with e(-) the usual Heaviside step 
function. 

Summarizing, to calculate the correlation and response we need P{E,t^), i.e. the 
solution of the original Barrat and Mezard model without a field, and M{t — t^\E^) 
from (fTT|) . The dependence on k is only through the latter. 
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4. Long-time scaling 

From now on we consider mostly an exponential density of trap energies, p{E) = exp{E) 
for E < 0; the glass transition temperature is then Tg = 1. We will also focus on the 
long-time scaling limit where memory of the initial conditions is lost and typical trap 
depths are large, \E\ ^ 1. The exit rate from such deep traps is 

r{E) = f dE'- = r'-r^ ^ = (19) 

^ ^ J-oo 1 + e/5(^'-^) Jo 1 + zl^ Jo 1 + z^ ^ ^ 

with c = 7iT/ sin(7rT). For Metropolis rates one has similarly 

rfE'e^'+ / ciE'e^'e-'^(^'-^)^ / dE' e""' + dE' e""' e'"^^'-^^ (20) 

-oo JE J—oo JE 

giving T{E) = c^e^ with 

= 1 + ^ = ^ = ^ (21) 
^ 13-1 1-T ^ ' 

If the timescales in the system are set by its age at long times, then typical values of 

T{E) - and therefore of - at time t should be of order t~^. This suggests the scaling 

ansatz P{E,t) = e^tV{e^t) where V{u}) is the normalized probability distribution of 

uj = e^t. The zero-field master equation Q for P{E, m, t) gives after integration over 

m 

d_ 
dt' 

Inserting the scaling form for P and neglecting the upper cutoff £" = in the integral 
transforms this into an integro-differential version for V{uj), 

u;V'{u;) = -(1 + c^)P(^) + / duj' ^ ^^^'^^.y (23) 

For Metropolis rates one only needs to replace c by and [1 + {uo/uj'YY^ by 
minKcu/cij')"^, 1}, and all statements for Glauber rates below can be translated to the 
Metropolis case in an analogous way unless specified otherwise. 

From it follows directly for — that P(0) = 1, as also found by Bertin |22] 
who used l/cj as the scaling variable. For large cj, on the other hand, i.e. relatively 
shallow traps, one expects effective equilibration and therefore V{uj) ~ e~^^ ~ uj~^ . 
(There is no density of states factor here because the density of states with respect to 
e is uniform on the allowed range < < 1.) At T = this power-law tail becomes 
an exponential as we will see below. For a numerical solution at T > 0, it is useful to 
rewrite ()23p as 

.Vi.) = f ..7(->— ■> /(-) = ji^' j^^, (24) 

For later we also note the following relation between successive moments of V{uj), 
obtained by multiplying ()23j) with cj" (—1 < n < j3 — 1) and integrating over u: 

fd^--»vM_j^_ .^^y'^^ (25) 



-P{E, t) = -T{E)P{E, t)+ J dE' p{E)w{E - E')P{E', t) (22) 



JduJuj''V{uj) V J l + zl^J c-c 
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2" vrT 
1 + zP sm [vTi [n + IjJ 

For Metropohs rates the analogous expression for the moment ratio ()25|1 is n/ (c'^ — cjf ) 
with 

cM = „i„{.-^l} = ^_^^|_^ (27) 

We next turn to the long-time behaviour of [j,{E, t — t^\E^). By arguments similar 
to those above, this should have the scaling form 

=7W(7,r), 7 = e^/e^-, r = (t - t„)e^" (28) 

where ^(7, r) is the normalized distribution of 7 after the rescaled time-interval r. The 
scaling form for M then follows directly as 

M{t-t^\E^) = JdE fx{E,t-t^\E^) = Mir), Mir) = J djU{-f,T) (29) 

The dependence on r = (t — t^)e^'" makes sense because, for \E^\ ^ 1, e^"' sets the 
scale of the exit rate from traps of depth E^. The dynamical equation (fTT|) becomes in 
the scaling regime 

M'ir) = -cMir) + nJ^-^^^ (30) 

This is similar in form to (|23|) for V{u!); accordingly, there is again an alternative version 
suitable for numerical iteration, 

e-" + ./V«Me-<-), «W^/^^^^ (31) 

The scaling of the correlation and response functions can now be deduced. For the 
correlation Q one gets 



r 



C{t,K)=C{{t-U)/U), C{x) = jdujM{xuj)V{uj) (32) 

which displays the expected simple aging scaling with x = {t — t^)/t^, implying that 
relaxation time scales grow linearly with the age. Similarly, the scaling form of the 
response function (fT^ is R{t, t^) = t~^lZ{x) with 

Tn{x) = {1-k) J dujMixuj)ritu) (33) 

= TcUoViuo) + [ duo' ; ^TT^ — } — , , (35) 

where we have used the integral 

duj' f dzz^ ^ f . f d\ 1 ^ . . 

^ 7T-T^ = T^ dzz [-—]-—-, = Tcu (36) 



[1 + {uj/uj'Y][1 + {uj'/ujY] J {l + zPf J \ dz)l + zP 
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ft 1 nX rjry.1 

X{t, U) = dt'- mt - t')/t') = x{x), x{x) = / T^7^(x') (41) 
Jt^ t' Jo 1 + x' 



For Metropolis rates, one has similarly 

r{uj) = J doo' [V{uj)e{uj' - uj){uj'/uj)-^ + V{uj')Q{uj - uj'){uo /uj')-^] (37) 

= Tc^ujV{uo) + J"duj'V{uj'){uj/uj')-^ (38) 

We will need the asymptotic behaviour of r{uj) below. For small u, one can approximate 
V{u!) in by 1. But the denominator ensures that a; and u' are of the same order, 
so the same argument can be applied to V{uj'). This gives 

r{uj) = 2Tcu (39) 

in the limit of small uj; for Metropolis rates one finds r{ij) = 2Tc^lj instead. For large 
uj, on the other hand, where V{uj) ~ one has 

giving the scaling r(a;) ~ uj^~^] the Metropolis case again differs only by c ^ c^. 
The susceptibility or step response x{t-,t^) = jl^dt' R{t,t') follows from (jHSj) as 

lo l + x' 

and again depends only on x = (t — t„)/t„, exhibiting simple aging scaling. Finally, the 
fluctuation-dissipation ratio (FDR) scales in the same manner, 

^ wiMu) ^ ^'^>- ^ - <^^' 

We recall that X = 1 corresponds to equilibrium, where the fluctuation-dissipation 
theorem holds; out of equilibrium, it then makes sense to define Teg = T/X as an 
effective temperature |TH1 dl 1201 1211 • This quantity allows a straightforward physical 
interpretation only when it is time- independent, at least within a given time-sector; 
the one of interest here is t — ~ t^. Since the definition (jl^ can be written as 
X = —{dTx/dt^)/{dC/dt^), such time-independence corresponds to a straight-line FD 
plot of x(t, t„) versus C(t,t^). More generally, if is used as the parameter varying 
along the curve [3 123 121| then X/T is the negative slope of such an FD plot. 

5. Standard Barrat and Mezard model at T > 

We first consider the case k = 0, i.e. the standard Barrat and Mezard model, at arbitrary 
temperature T (< Tg = 1) and with an exponential density of states. One then 
has, from ^ and M{t) = exp(-cr). The scahng forms ^ and ^ 

of correlation and response function thus simplify to 

C{x) = Jduje-'''''V{uj) (43) 

and 

Tn{x) = fdue-^'^'^r^uj) (44) 
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From this the large-x asymptotics can be obtained directly: only the region ~ l/x 
contributes to the integrals so that in C(x) we can set V{uj) ~ 'P(O) = 1, giving 

C(x) = — (45) 

cx 

as found previously by Bertin 22j. For TZ{x), using leads similarly to 

2T 

Tn{x) = — (46) 

The FDR from ()42|) thus behaves asymptotically as X{x) = 2T/x. In terms of C this 
gives X = —Tdx/dC = 2TcC for small C. The limiting FD plot therefore always 
starts parabolically at the top, x{C) = x(0) — cC^ for C ^ 0. The FDR itself decays 
to zero asymptotically, X{x oo) = for all T < Tg, as in the case of the Bouchaud 
trap model IH IHI- 

More interesting is the short-time FDR X{x 0). Naively, this is from (j^^ 

^ ' -C'(0) JdujcuViu) cJduuViu) ^ ' 

using (jini) and with x = as well as and (or for the Metropolis 
case). The effective temperature associated with this FDR is therefore = 1/2, 
independently of T. However, this conclusion only holds for T < 1/2. For larger T, the 
integral 

- C\x) = J rfw c we-'=^'"P(cj) (48) 

is actually divergent for x = because of the uj~^ tail of V{uj). For small but nonzero 
X, the exponential acts as a cutoff at ~ 1/x so that — C'(x) ~ (l/x)^~^ = 
in agreement with the short-time singularity in C identified in [22] ■ The response 
TTZ{x) = Jduj e~'^^'^r{uj) has the same small-x singularity because r{u)) ~ cj^"^ for 
uj oo. In fact, for small x one can replace r{uj) by its asymptotic form ()4()j) to get 



T7^(x) = J ducue-^'^'^Piuj) = -C'{x) (49) 

This shows that the short-time FDR is, for T > 1/2, 

X{x ^ 0) = 1 (50) 

In this temperature-range the FD relation between x cind C is therefore always of a 
pseudo-equilibrium form in its initial part, with a (negative) slope equal to 1/T. This 
matches continuously with the constant slope of l/Tes = 2 found above for T < 1/2. 

Intuitively, one can understand the occurrence of pseudo-equilibrium behaviour for 
T > 1/2 by looking at the average hopping rate. This is 

r(t) = jdET{E)P{E,t) = t-^ jdujcujV{uj) (51) 

where the second form applies in the scaling regime. For T < 1/2, the integral converges 
and r(t) ~ 1/t as one would naively expect from our scaling assumption: typical 
relaxation times are ~ t, thus typical rates are ~ 1/t. For T > 1/2, however, the 
integral is divergent at the upper end and one has to take into account the cutoff at 




Figure 1. Numerically calculated scaling distributions V{lo) for Glauber dynamics 
at T = 0.2, 0.4, 0.6, 0.8 (bottom to top on the right). The dashed lines indicate the 
expected asymptotic scaling 'P{ijj) ^ lo^^ . 

E = 0, corresponding to u = t, leading to r(t) ~ t^~^. This is entirely dominated by the 
very small probability of being in shallow traps with atypically fast relaxation rates of 
0{1). Now the response R{t,t^) and the initial decay {dC/dt^){t,t^) of the correlation 
function are sensitive only to hops taking place between and t. For small t — 
these are precisely the same events that dominate T{t). Since they are in the "effective 
equilibrium" tail of V{uj), it is then not surprising to find a pseudo-equilibrium form of 
the FD relation. 

Fig. ^ shows some numerically calculated scaling distributions V{uj) for Glauber 
rates; these exhibit the expected asymptotic behaviour ~ for large u. To avoid 
end effects in the iterative numerical solution of (j2H), we stored not V and / themselves 
but V{uj){l + u)^ and similarly /(a;)(l + These functions have nonzero limits 

f or LJ ^ oo (and a; — >^ 0) and so are suitable for evaluating the required integrals over 
a; = . . . cxD. 

Fig. 121 displays the resulting FD plots of x vs C; these are valid in the limit of long 
times which we have already taken by working in the scaling regime. (To get reliable 
results for x(x), one has to do the x'-integration from ()41|) before the tu-integration 
in (jnSl); otherwise the singularity of 7l{x) for a; — > at T > 1/2 leads to problems.) The 
initial slopes agree well with the theoretical predictions, as shown by the dashed lines. 
The asymptotic slopes for C ^ are likewise consistent with the predicted value of 0. 
The inset of the figure explores this region in more detail, showing (1 — x)/(cC^). From 
the analysis above this quantity should converge to 1 for C ^ 0. The numerical data 
are consistent with this, though for T = 0.8 the approach to the limit is very slow. This 
makes sense: as T tends to Tg = 1 from below, the FD plot approaches a straight line. 
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Figure 2. Numerically calculated long-time FD relations x(C') in the original Barrat 
and Mezard model (k — 0, Glauber dynamics). Solid lines are for T = (exact 
0.2, 0.4, 0.6 and 0.8 from top to bottom; the results for T = and T = 0.2 are almost 
indistinguishable. The dashed lines show the theoretically predicted initial slopes (2 
for T < 1/2, and 1/T for T > 1/2). The inset graphs (1 - x)/{cC'^) versus C; our 
theory predicts that this quantity converges to 1 for C ^ 0. 



and so the quadratic expansion around C = will be valid in a region that shrinks to 
zero in the limit. 

One notable feature of the numerical results is that the asymptotic value x{C = 0) 
is equal to 1 for all T < 1, at least to within our numerical accuracy of around 10^^. 
This is as in the Bouchaud model JU]; and has been conjectured also for the Barrat and 
Mezard model, on the same physical grounds [TT]: the susceptibihty freezes to its value 
at T = Tg = 1 as T is decreased below the glass transition. It ought to be possible to 
confirm this result analytically, but we have not yet found a way of doing this. 

It should be noted that the FD plot is not quite T- independent for T < 1/2. While 
such T- independence had been suggested by the simulations of [T7], the scaling- regime 
numerics shown in Fig. El clearly rule it out. Also, after a little reflection one sees that 
higher derivatives of C{x) diverge for x ^ already at lower temperatures; e.g. C"{x) 
diverges for T > 1/3, C'"{x) for T > 1/4 etc. 

6. Slowly decorrelating observables, T = 

In this section and the next we consider slowly decorrelating observables, i.e. k > 0. 
We begin by analysing the zero temperature dynamics, where a number of results can 
be derived for a general density of states p{E) and finite times (i.e. without taking the 
scaling limit). 



Trap models with slowly decorrelating observables 



12 



Compared to the standard 0, the additional task is to calculate 

lj{E,t-t„\E^). At T = 0, one has w{E - E') = Q{E' - E) and V{E) = jf^dE' p{E') 
for both Glauber and Metropolis rates; equation ()14|1 thus simplifies to 
s + T(E) 1 /"O 

—-^-lf^{E,s\E^) - ——6{E - E^) = K / dE'fi{E\s\E^) (52) 
p{E) p{E„) Je 

Differentiating w.r.t. E gives a differential equation for fi{E,s\E^) as a function of E. 
This is easily solved, with the result 



jj,{E,s\E^ 



d 



s + r(E„) dE 

1 d 



Q{E^-E)exp l^K 

s + T{E. 



dE' - 

E S 



P{E') 



V{E') 



QiE^-E) 



T{E) 



(53) 



s + T{E^) dE 

Before we exploit this to obtain M and thence the correlation and response functions, it 
is worth noting from ((7j) that, for the limiting case k = 1, fi{E, t — t^\E„) is simply the 
propagator P{E,t — t^\E^) of the original Barrat and Mezard model without a field. 
Inverting the LT in thus yields the exact zero-temperature propagator 



P{E, t - U\E„) = \e{E^ - E)e 

oE L 



-V{E){t-t^) 



(54) 



Applying this to 

P{E,0) 



gives as the general solution starting from an initial distribution 



CdE'P{E\ 0) 

JE 



(55) 



At long times most of the mass of this is a low E, where the integral inside the 
square brackets can be set to one. One thus recovers P{E,t) = p{E)texp[—T{E)t] 
as the long-time scaling form of the distribution over trap energies, independently 
of the initial distribution For an exponential density of states, p{E) = , this 
becomes P{E,t) = e^tV{e^t) with the scaling function V{uj) = exp(— cu). This has an 
exponential tail as anticipated above, and one easily checks that it solves the T ^ 
limit of (121 • 

Returning now to the result (fS^ for general k, we obtain for the LT of M(t — t^\E^) 

1 f s + T{E„) y 
' s + T{E^) [ s ) 

Inverting the LT gives a hypergeometric function. 



M{s\E^) = jdEfi{E,s\E^ 



s-^'is + ViE^ 



IK-l 



(56) 



M{t - K\E^) = iFi 1 - /€, 1; -T{E^){t - U 



(57) 



If we now specialize to the exponential distribution of trap energies, we have T{E) = 
and so (j57jl is exactly of the scaling form (j29|l discussed above, with 



M(r) = iFi{l-K,l;-T 



E 

fc=0 



-r 



{51 



(Non-integer factorials are defined in terms of the Gamma function, a! = r(a + 1); we 
avoid the use of F for such constants to prevent confusion with our exit rates T{E).) 
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Figure 3. Long-time FD plots x(C) at T = 0, for k ^ 0, 0.2, 0.4, 0.6, 0.8, 1 from 
bottom to top. 

Both of these results can also be obtained directly from the relevant equations (fTT|l 
and (pn)) for M(t\E„) and A4{t), by representing these quantities as power series in t 
and r, respectively, and determining the coefficients recursively. 

From A^(r) we get the scaling limit of the correlation function, using V{u!) = 
exp{-uj), 

C{x) = / dijM{xij)V{u;) = V ) .. 'A -xf = (1 + xf-^ (59) 

so that C(t,t^) has the simple long-time form C{t,t^) = {t^jt)^^'^. This result can also 
be derived without recourse to series expansions: denoting M.{-) the LT of M.{j), we 
have from ()32|) 

Jdxe-''''C{x) = JduJUJ-^M{a/u)e-'' = j dxx-^M{x~^)e-'''' (60) 

so that C(x) = x~^M{x'^). With M{x~^) = (a;"^)~''(x"^ + 1)"^^ from dSHI) one gets 
C{x) = (1 + x)""^ as before. 

For the response, we note that the denominator factor in (jH^ tends to Tlj6{lj — u') 
for T — > 0; the same conclusion also holds in the Metropolis case pTjl . This gives for 
the scaling function 

7^(x) = 2(1 - /t) J duo M{xuj)ujV{uj) = 2(1 - ^^^^^yj^i i'^)'' (61) 

= 2(1 - k)-^[xC{x)] = 2k{1 + x^-^ + 2(1 - n^il + x^-^ (62) 
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For the susceptibility this imphes, from (j41|) . 

Xix) =2k\1-{1 + xr~'] + ^IIZ^ [l - (1 + xr-^] (63) 

= 2f,[l-C] + ^ilz^ [i _ (64) 

where in the second equahty x has been expressed in terms of C to give the long- 
time FD relation x{C)- This is displayed in Fig. 01 for some representative values 
of K. The slope of the FD plot is -x'{C) = 2k + 2(1 - k)C1/(i-«), with initial 
value — x'(C* = 1) = 2 independently of n, corresponding to an "effective short-time 
temperature" of T^s = 1/2 (= Tg/2 in dimensional units). The asymptotic slope, on 
the other hand, is — x'(C = 0) = 2k and depends continuously on k. In the limit k — 1 
(where C{x) = (1 + x)'^^^ decays very slowly with x, and x{x) grows correspondingly 
slowly), the FD plot becomes a straight line of negative slope 2, suggesting that this 
slowly decorrelating observable measures a well-defined T^s = 1/2. This result lends 
support to the straight-line FD plots found in the simulations of Junier and Bertin jl7j . 

For later, we note that by expressing the ratio of factorials in as a Beta function 
integral and then performing the sum over k one gets an alternative form for A1(r), 

Mir) = -T^ — [^dze-'^z-^{l - zY~^ (65) 



. -k)!(k- 1, 

This implies in particular that A^(r) = r^'^^^'^y^K, — 1)! for large r and k, > 0. 
7. Slowly decorrelating observables, T > 

Finally we consider the most general case of dynamics at nonzero temperature and 
observables with incomplete decorrelation (k > 0). We focus directly on the scaling 
limit for a system with an exponential density of states; the key issue is again solving 
for A^(r). For short times this is easy: from (j3Up . A^'(O) = —c + kc and thus 
M{r) = 1 - (1 - ft)cr + C(r^). From this the initial FDR for T < 1/2 can be worked 
out as in (|T7|) 

^ ' -C'(0) jdu{l- k,)cujV{uj) ^ ' 

The K-dependent factors (1 — /t) cancel, so the initial slope of the FD plot remains 2 for 
all K. Metropolis rather than Glauber rates again only replace c by everywhere. 

For T > 1/2, the linearization in r breaks down as before, and the integrals for TZ{x) 
and 1 — C{x) will be dominated by large ~ 1/x for x 0. One can therefore use the 
asymptotic behaviour of V{uj) oc uj~^ and r{uj) oc cuj^~^ (with the same proportionality 
constant, according to (jlOI)) to write 

1 - C(x) oc Jduj[l-M {xuj)]uj~(^ = x^-^ jdT[l- M (r)]r-^ (67) 
T7^(x) oc{l-K)cJduM{xu)u^-'^ = {l-K)cx^-^JdTM{Ty-^ (68) 
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These singularities (but not their prefactors) are the same as for k = 0, and 
correspondingly the FDR again approaches a nonzero limit for x 0: 

{l + x)C'{x)^{(3-l)Jdr[l-M{T)]T~(' -JdTM'iry-f' ^ ' 

At first sight the value of X{x — > 0) appears to depend on the precise functional form 
of M.{t). But in fact, by multiplying (IHUj) by r^~'^ and integrating over r one deduces 

J drM'iry-^ = -cJdTM {ry-^ + k j dr' M (r') j ^ ^ ^[^^i^y C^O) 

The last r-integral evaluates to c(r')^~^ and this implies from ()fi9j) that A'(x 0) = 1 
for all K (and T > 1/2). Physically, this result — which can be derived similarly for 
Metropolis dynamics — is supported by the same intuition as for k = 0: the initial 
response and decay of the correlation, and hence the FDR, are dominated by the very 
small fraction of histories which pass around time through the shallow traps near 
the top of the energy landscape {E = C(l)), where an effective pseudo-equilibrium is 
established. 

What about long time intervals, x ^ 1? From the definition of the dynamics, it 
is clear that C{t,t„) is the average of = e^^^'^ over the distribution of the number 
of hops j between t„ and t. For k ^ 1, a large number j ~ l/(— In/t) ^ 1/(1 — k) 
of hops is needed to get any significant decorrelation. Across the corresponding long 
time intervals the number of hops should average out, leading to the naive prediction 
C ~ exp[(lnfi:) J^^dt' T{t')]. Using (jSH), this gives the estimate 

C{x) = {1 + xy^^-'^l r] = cJdujujV{u) (71) 

Intuitively, one expects that similarly A4{t) ~ r^'''ii'^-'^) for large r. The constant r] can 
be evaluated explicitly, by taking the n — > limit of (|^. giving 

^ ^ ^ sin^(7rT) ^ tan(7rT) ^^^^ 
(ttT)^ cos(7rT) ttT 

The analogous expression for the Metropolis case is, using (fTTj) . 

T] =c I duj ojV ioj) = — — T = — = — 73 

' J ^ ' (3-l(3{(3-2) (3-2 1-2T ^ ' 

and together with (ffT| confirms the estimate for the correlation function scaling in 
eq. (22) of ^7j. Both t] and r]^ diverge as T — 1/2 from below, signalling a breakdown 
of the above naive reasoning for higher temperatures; we return to this point below. For 
T —* 0, 1] and r]^ both approach unity so that M.{t) is predicted to decay as t~^^~'^\ 
exactly as we found by exact calculation for T = 0. 

To find the asymptotics of A^(r) for general T and k, one substitutes the ansatz 
Ai{T) ~ r"'*' into (jHUj) . The l.h.s. is then subleading and the leading terms on the r.h.s. 
have to cancel, giving the condition 

C = KC_A (74) 
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Figure 4. Calculated values of the decay exponent A of A^(t) ~ t as a function of 



K. The curves are for Glauber dynamics and T 
top. 



0, 0.2, 0.4, 0.6, 0.8 from bottom to 



for A, or equivalently fi;sin(7rT) = sin[7rT(l — A)]. For c_a to be finite one needs A < 1, 
and A^(r) decreases with r so A > 0. In this range of A the condition ()74j) has only one 
solution, 

arcsin[K sin(7rr)] 



A 



vrT 



(75) 



For small k this exponent approaches 1, with A = 1 — k/c + see Fig. |3 The 

behaviour for k 1, on the other hand, depends strongly on temperature. For T < 1/2, 
A = (1 — k) tan(7rT)/(7rT) = rj{l — K) to first order in 1 — k, exactly as the naive argument 
had predicted. For higher temperatures T > 1/2, the exponent approaches a nonzero 
value A = 2 — /? in the same limit. This is rather striking since at k = 1 directly one 
expects A^(r) = 1 and hence A = 0: the limit k ^ 1 is discontinuous for T > 1/2. We 
discuss the physical reasons for this at first sight surprising behaviour below. 

The above predictions for A apply to the Glauber case. Using the same arguments 



for Metropolis rates, one is led to the condition 

f3/{/3 - 1) = «:/5/[(l - A)(/5 + A - 1)], giving 



KC 



M 



A 




<P - 1) 



This reads explicitly 



(76) 



Similarly to the Glauber case, A ^ for K^latT<l/2, with leading order behaviour 
A = (1 — K){j3 — l)/(/3 — 2) in agreement with (fTSj) . For T > 1/2, on the other hand, 
A 2-/9 and this limit value is the same as for Glauber dynamics. In general, however, 
as soon as the decorrelation of the observable in a jump is not perfect {k > 0), the decay 
exponents (fTSj) for Glauber and (fTBj) for Metropolis rates are different. This provides 
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an interesting example where microscopic details of the assumed transition rates affect 
not only quantitative details (prefactors) but also qualitative aspects (exponents) of the 
dynamics. 

The power-law decay of A^(r) implies similar large-x asymptotics for the response 
and correlation functions, consistent with the naive expectation explained above: 

C{x) ocx-^ Jdujuj-^V{uj) (77) 

Tn{x) oc (1 - k)x'^ J duj uj-\{uj) (78) 

In particular, for T > 1/2 and k — > 1, the correlation decays as C{x) ~ x^'^ with 
A = 2 — /3. This exact value is consistent with the argument by Junier and Bertin |^ 
that the decay exponent A should be close to T = 1//? for T ^ 1; to first order in 1 — T 
this estimate agrees with the exact prediction. It deviates progressively for lower T, 
however, with the exact value \ = 2 — (3 dropping to zero at T = 1/2. This is as it must 
be by continuity, since for T < 1/2 we found that A vanishes for k 1. As a check of 
our theory we have also compared the predicted exponent (fTBjl for Metropolis dynamics 
with the numerical correlation function results from Fig. 9 of |T7j for k = 0.9, and found 
very good agreement. 

From the above large-x behaviour of C and TZ we can deduce the asymptotic 
FDR dig) as 

xi. - oo) = '^-;>^"Xrt' (™> 

Using one finds that the numerator integral equals T[c+ (1 — A)c_a] / du uj^^^V{uj). 
The remaining moment ratio can then be calculated from with n = —A to give 

_^ (i-.)r|c + (i-A)c.>| 

C-A - C 

Inserting the condition (f7^ determining A simplifies this to 

X{x ^ oo) = T{1- \ + k) (81) 

One can check that this expression is valid also for Metropolis rates, with the value of A 
then given by (fTUj) instead of A. For both Glauber and Metropolis rates the asymptotic 
FDR X{x —>■ oo) approaches zero for k ^ 0, consistent with the results for the ordinary 
Barrat and Mezard model from Sec. |S| The behaviour for n —>■ 1 again depends on T. 
For T < 1/2, A and thus X{x oo) ^ 2T. For T > 1/2, on the other hand, 
A 2-/9 and X{x oo) — > 1. Both of these values are identical to the (K-independent) 
short-time FDR X{x 0) found above. This strongly suggests that the FD plots are 
straight lines for k, ^ 1, both above and below T = 1/2. We now proceed to show this, 
by considering the FDR X{x) for general x. 

To this end, it will be useful to have an expression for C'{x) that is similar in form 
to 7^(x). Starting from (jSl one has C'{x) = Jdu M'{xuj)ujV{uj). With (jUl) this can be 
written in terms of itself as 

C'{x) = -c fduM {xu)uV{lo) + k fduM {xu) f du' (82) 
J J J 1 -|- [u) j uj y 
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The a;'-integral is, from (|23p . {d/duj)[ujV{uj)] + cujV{uj). The second term in this sum 
gives a contribution of the same form as the first term on the Lh.s. of (jH^ . and after an 
integration by parts 

C'{x) = -c(l - 1^) J du M{xuj)ujV{uj) - kx j du M'{xuj)ujV{uj) (83) 

The last term is now just —k,xC'{x) and we end up with the relatively simple expression 

- (1 + Kx)C'{x) = c{1-k) JdcuM {xuj)ujV{uj) (84) 

which also holds for Metropolis dynamics (c c^). The FDR ()42j) then becomes a 
ratio of integrals involving only A^(-), rather than A4{-) and as before, 

y(r-] = 1 + Jduj Mixuj)r{uj) 

^ ' l + x cSduoM{xu)ujV{u) ^ ' 

We can now see how this simplifies for k — 1. For T < 1/2, 7W(r) — 1 as k — > 1 for 
any fixed r. The u;-integrals remain finite in this limit, so that 

^(.) - = 2T (86) 

c J dujujy[uj) 

exactly as in ()47|): A:'(x) = 2T becomes independent of x, and the FD plot is therefore a 
straight line of slope X/T = 2. For T > 1/2, on the other hand, A^(r) is asymptotically 
a power law r~'^, with A approaching 2 — (3 from above as k — > 1; thus for any fixed 
X > also M.{xuj) ~ uj~^ for large cj. But r{uj) ~ ti;'P(ti;) ~ and so the integrals 
in (|HH|l become divergent at the upper end as k — 1. They are therefore dominated by 
large cj- values, where r{uo) = cujV{uj). This shows that X{x) 1: the FD plot is again 
a straight line, but now of equilibrium slope 1/T. 

With similar arguments we can also work out the correlation functions for k — > 1. 
Consider T < 1/2 first, and divide (j5^ by C{x) to get 

n , _,C'{x) _ c{l-n)^dujM{xuj)ujV{uj) 

^C{x) ^dujM{xuj)V{uj) ^ ^ 

In the limit k ^ 1 we can set ^ 1 again and get — c?(lnC)/(i(ln(l + x)) = 
c(l — k) jdujujV{uj) = r](l — k). This implies C{x) = (1 + x)~^^^~'^\ exactly as the 
naive argument (fTTj) suggested 

For T > 1/2, consider the integral on the r.h.s. of ()84|1 . This becomes divergent at 
the upper end for /t — > 1 as explained above, and dominated by large u. J^{xu) can 
therefore be replaced by its asymptotic form ~ (xcu)"'^, and the r.h.s. of becomes 
proportional to which in the limit k — 1 is (The remaining integral diverges 

as K — > 1, but combines with the 1 — k prefactor to give a finite limit.) Thus, for T > 1/2 
and K — > 1, 

. a;^"' ^ sin[7r(/? - 1)1 /- , , (x'^"^ 

- C'(x) oc ; , C(x) = — ^ / dx' \ ' , 88 

In the expression for C(x) we have explicitly put in the proportionality constant, which 
follows from the equal-time value C(0) = 1. The same result also holds for Metropolis 
dynamics, since both the limit value of A for k — 1 and the asymptotic behaviour of 
V{u}) and r{uj), on which the argument relies, are the same. 




Figure 5. Calculated Mir) for T = 0.4 and k = (exact), 0.2, 0.4, 0.6, 0.8, 0.9, 0.95, 
from bottom to top. Dotted lines indicate the predicted r"'^ asymptotes. 

7.1. Numerical results 

We solved numerically to get 7M(r), using a representation which takes into account 
the asymptotic behaviour t~^, and then evaluated correlation and response. The results 
for A4{t) at T = 0.4 are shown in Fig. and do approach the predicted asymptotes. In 
particular, the slope A of the asymptotic power law decreases to zero for k ^ 1, while 
for T = 0.6 (Fig. ^ it approaches the nonzero limit 2 — (3. The dashed line shows the 
limiting form of M.{t) for k — > 1, which for T > 1/2 is nontrivial. Note that we focus on 
Glauber dynamics throughout this section; graphs for the Metropolis case would look 
broadly similar as discussed above. 

Fig. [7| shows the scaling function C{x) of the two-time correlation at T = 0.4, again 
for a range of values of k. As expected, the asymptotic decay follows the same power 
law as for A^(r), i.e. C{x) ~ x~'^. In the inset we show C(x)/(1 + x)~'^ to demonstrate 
that this approaches unity in the limit k ^ 1. Fig. |H1 displays the analogous data 
for T = 0.6. The predicted asymptotic behaviour is again observed, but now C{x) 
approaches the nontrivial limiting form for k, —>■ 1 predicted by (jHHj) (dashed line). 

The resulting FD plots (Fig. IHl left) for T = 0.4 move upwards with increasing k 
as for T = and are consistent with the approach to the predicted straight line of slope 
2. For T = 0.6 (Fig. El right) the overall trend is similar, but the limiting straight line 
now has slope 1/T, again as predicted. 



Trap models with slowly decorrelating observables 



20 




0.01 



0.001 



Figure 6. Analogue of Fig. [S] for T = 0.6, for the same values of k. The dashed line 
shows in addition the numerically calculated limiting form of Ai{T) for k ^ 1. 




0.01 



Figure 7. Correlation function C{x) for T = 0.4 and the same k as in Fig. |S1 Dotted 
lines indicate the predicted asymptotes x^^. The inset shows C(x)/(1 + x)^^ versus 




Figure 8. Correlation function C{x) for T 0.6 and the same k as in Fig. (SJ the 
dashed line shows the predicted k ^ 1 limit (|88|l . Dotted lines indicate the predicted 
asymptotes x^'^. 




Figure 9. FD plots for k = 0, 0.2, 0.4, 0.6, 0.8, 0.9, 0.95 1.0 (the last one being exact) 
from bottom to top, at temperature T — 0.4 (left) and T = 0.6 (right). 



7.2. Discussion 

The results found above for the hmit k 1 of slowly decorrelating observables provide 
strong support for the arguments of Junier and Bertin |T7j , in particular in terms of the 
straight-line FD plots and the asympotic decay of C{x) for T < 1/2. In addition they 
clarify the limiting behaviour of C{x) for T > 1/2, both in terms of the asymptotics 
(~ x^""^) and the full (and rather simple) functional form ()88|1 . Finally, our theoretical 
analysis predicts how the decay exponent A interpolates between its value for k, = 
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(standard Barrat and Mezard model) and the k = 1 limit of slowly decorrelating 
observables, and reveals that the details of this do depend on the choice of microscopic 
transition rates (Glauber versus Metropolis). 

One of the most surprising findings is the discontinuity in the approach to k = 1 for 
T > 1/2: strictly at k = 1 the correlation function cannot decay and so C{x) = 1 and 
M.{t) = 1. For K just below 1, on the other hand, we found that these scahng functions 
do decay for finite values of their arguments, and approach nontrivial limit forms for 
K — > 1. At first sight this seems strange: as explained above, C{t,t^) is the average 
of over the distribution Pj{t,ty^) of the number of jumps between and t, and as 
K ^ 1 this average must approach unity. Indeed, this statement is true for any fixed 
pair of times t and t^. The discontinuity only arises in the limit of large times, where 
Pj consists of two essentially separate parts. The first contains finite j, corresponding 
to histories where the trap energies remain of the same order as the initial one: the 
system wanders among deep traps. The contribution to the correlation function arising 
from this part of Pj is indeed continuous for k — > 1. However, Pj also has a second part 
(containing a finite fraction of its probability mass) at values of j which are divergent 
with tyf, presumably - for consistency with the scaling of the hopping rate - as t^^. 
This part arises from histories which make an excursion from deep traps to the top of 
the landscape. In the limit t^ ^ oo, the correlation decays to zero for any such history, 
as long as K < 1 (i.e. for any small but nonzero 1 — /t). This is why the scaling functions, 
which arc calculated in precisely this long-time limit, are discontinuous at k = 1. For 
finite t„ one expects to sec the crossover between correlation functions which essentially 
do not decay, and the limiting forms calculated for k — 1, to occur for values with 
K*^ of order one, i.e. 1 — k ~ t^~'^ for large t^. 

This argument also tells us how to interpret the k — > 1 hmit of C{x): it is the 
probability of never having escaped to the shallow traps at the top of the energy 
landscape. It is rather intriguing to see that this can be calculated exactly while the 
K — > correlation, i.e. the probability of not having jumped at all, cannot as far as we 
are aware. 

Finally, it is interesting to see what happens to our naive estimate from above, 
C ^ exp[(lnK)//^dt'r(t')], for T > 1/2. This would give C(t,t„) = exp[-const x (1 - 
n){t^~f^ — t'^f^)] which decays on subagiug timescales t — t^, ^ ^w^^- Applied to the 
entire correlation function this is plainly wrong: the correlation function must decay 
more slowly as k increases, but ior k — we have simple aging, so k > must also 
give simple aging or an even slower decay. We now see that the argument only applies 
to the decorrelation caused by large (diverging with t„) numbers of jumps, where it 
predicts correctly that the resulting contribution to the correlation function vanishes 
in the scaling limit. The remainder of the correlation function remains finite (and has 
simple aging), however, because it relates to the jumps among the deep traps which 
happen at rate ~ 1/t, rather than to the total number of jumps which is dominated by 
the shallow traps. 
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8. Bouchaud model 

The dynamics of the Barrat and Mezard model can be interpreted as being slowed 
down by entropic barriers at low temperature: the Glauber (or Metropolis) transition 
rates penalize large energy increases, so that the system is forced to search for lower 
energy states of which there are an ever decreasing number. To complete our analysis, 
we now compare with the Bouchaud trap model jl], where jumps take place with rate 
exp{(3E) independently of the energy of the arrival trap; thus almost all jumps return 
the system to the top of the energy landscape. Slow dynamics then arises purely from 
activation effects, i.e. the decreasing rate with which jumps from increasingly deep traps 
can take place. We will be interested to know in particular whether slowly decorrelating 
observables again produce straight-line FD plots. 

Applying the intuitive reasoning from the previous section to the Bouchaud model, 
we see that once the system has returned to the top of the landscape it will make a 
large number of fast jumps, which will decorrelate the values of the observable m for 
any k, < 1. Thus, the correlation function C{x) should be completely independent of k 
for < K < 1 since, in the scaling limit, correlations are maintained only if no jumps at 
all take place (j = 0). 

If the Bouchaud model, the natural prescription for the coupling to the field is to 
modify the escape rate from a trap from exp(/?£') to exp[P{E — hm)] IHIH]- This is just 
the original rate multiplied by exp(— /3/im), and one is led more generally to consider 
an arbitrary trap model with multiplicative rates [Ij ITT|: 

The notation here should be self-explanatory; the arrow points from the departure to 
the arrival trap. In the absence of a field all governing equations are as before, except 
for the replacement of w{E' — E) by w{E' ^ E). The master equation then reads 

^P(E, t) = -T(E)P(E, t)+ [ dE' p(E)w(E ^ E')P(E', m', t) (90) 
ot J 

with T[E) = JdE' p{E')w{E' ^ E), while the magnetization decay function obeys 

—M{t\E) = -T{E)M{t\E) + KjdE' M{t\E')p{E')w{E' ^ E) (91) 

and the two-time correlation function is given by Q as before. To get the response, 
one can start again from ()15|1 . The change in P{E,m,t^ + At) from its value without 
a field is given by the direct analogue of (fT?)|) . 

AP(E, m, t„ + At) = At y" dE^ dm^ x 

[p{E)p{m\m^)Aw{E,m ^ E„,m^)p{m„)P{E^,t^) 

-p{E^)p{m^\m)Aw{E^,m„ <— E , m) p{m) P {E , t^)] (92) 

Now to linear order in h, Aw(i?,m ^ E„,m^) = Ph[{l — Qm — (m^]w{E 4— E^). 
Inserting this and integrating over m^, and then over m in p5|) . yields 

TR{t, t„) = JdE dE^ M{t - t„\E) {{-(k + 1 - C)p{E)w{E ^ E^)P{E„, t^) 



Trap models with slowly decorrelating observables 



24 



+ [C - (1 - Ok]p{E^)w{E„ ^ E)P{E,t^)} (93) 

As in the case k = one can now relate the response function to derivatives of the 
correlation function. From (jH)) together with ()91|) one has 



d_ 

dt 



C{t,t^) = JdE 



T{E)M{t-t^\E) 



+ K JdE' M{t -t„\E')p{E')w{E' ^ E) 



P{E,t^) (94) 



For the derivative w.r.t. the earlier time t^ one gets similarly, by also using (j9(Jj) . 

—C{t,t^) = {l-K) fdEdE' M{t~t„\E')p{E')w{E' ^ E)P{E,t^) (95) 

This is proportional to the first term in the expression (j93|l for the response, while the 
second term in (j93p is proportional to the first one in (j94j) . bearing in mind the definition 
of T{E). Altogether one gets the simple relation 

TR{t,t^) = (1 - 0(1 + n)—Cit,t^) - [C - (1 - C)'^]g-^C{t,t^) (96) 

which reduces to the known result j3J for k = as it should and applies to all 

trap models with multiplicatively perturbed rates (jHU)) . At high temperatures, where a 
time-translation invariant equilibrium state is reached, dC/dt^ = —dC/dt. As required 
for consistency with equilibrium FDT, the coefficients of these two quantities in ()96p 
add up to unity. Out of equilibrium at low temperatures, an equilibrium FD relation 
is still recovered when the second coefficient vanishes, i.e. when C, = + k). It is 
easy to see that this is exactly the case where the exit rate T{E,m) from a given state 
is unperturbed by the field h (to linear order), just as for k = piU lllj. This implies 
that for slowly decorrelating observables (k 1) one only gets a straight-line FD plot 
when C = 1/2, i.e. when the hopping rates depend on the difference m' — m as they do 
for Glauber dynamics. 

For the Bouchaud model specifically, where w{E' ^ E) = exp{PE), equation 
for M(t\E) is easily solved by LT to give 

M{s\E) = ( 1 + , Gis) = fdE ^ (97) 

^ ' ^ s + ef'^ \ + KsGis) J ^ ' J s + ef'^ ^ ^ 

The scaling limit for P{E,t) is reached for large t and low E with te^^ of 0{1). This 
corresponds to small s, of order e^^. Since G{s) ~ s^~^ for small s, one sees that in this 
regime the second term in M{s\E) always becomes negligible compared to the first, as 
long as K < 1. This means that the K-dependence drops out, and M{t\E) = exp{—te^^) 
in the scaling regime. The correlation function only picks up these scaling contributions 
and is therefore ^-independent as expected. 

The same argument does not apply to the response, which from ()96|) is obviously 
dependent on n. The reason is clear from (j^Hj) : the contribution from M{t — t^\E) 
is weighted with extra factors of either p{E) or w{E^ ^ E) (= e^'^ in the Bouchaud 
model). This means that the behaviour of M{t — twl-^') for E = 0{1) is dominant, and 
this does depend on k. 
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8.1. Arguments for straight-line FD plots 

The Bouchaud trap model provides a useful point of reference from which to revisit the 
reasoning put forward in jT2j for the emergence of straight-line FD plots in the limit 
K, 1. This has at its core the idea of considering first the change of the magnetization 
m during j jumps, and then to average over the distribution of the number of jumps 
between t^ and t. There is no problem with this for the case of the correlation function, 
since in the absence of a field the magnetization is just being "convected along" with 
the usual trap model dynamics, without itself affecting the dynamics. 

For the susceptibility, on the other hand, the situation is somewhat more subtle. 
Consider a field switched to some nonzero value (and held there) at time t^,- Then one 
can write generally 

Q ^ r 

X = -gj^^ J dmdm„dE„mP{m\j,h,E^,m„,t -t„) x 
j 

X Pj{h, E^, m„, t - t„)P{E^, m„, t^) (98) 

where P{m\j,t — t^, h, E^,m^) is the distribution of m at the end of a time interval 
t — t^ of constant field strength h, given that j jumps have taken place during this time 
and the system started from a trap with energy E^ and magnetization m„. Pj{. . .) is 
the distribution of the number of jumps, which depends on the same variables. Finally, 
P{E^,m„,t^) is the distribution of trap energies and magnetizations at t„, which can 
be written as P{E^,t^)p{m^) if as we assume the field has been off up until time t^. 
Writing the m-average in as fhj{h, E„, m„, t — t^) one thus has 
d . f 

^ = ^ 51 _/ dm^ rhj{h, E^, m^, t - t^)Pj{h, E^, m^, t - t^)p{m^)P{E^, t^) (99) 

To complete the argument of Junier and Bertin fTP^ two assumptions now need to be 
made: (A) Tfij{h, E„, m^, t — t„) is the same as the average of m directly after the j-th 
jump starting from a state with E^ and m^. (B) The distribution of the number of 
jumps Pj{- ■ ■) is independent of m^, so that rhj{. . .) can be averaged separately over 
m^. The first assumption is in general invalid because the fact that there has not been 
a further jump between the j-th jump and time t favours magnetizations m of the same 
sign as the field h. Indeed, in the standard = 1, k = 0) Bouchaud model assumption 
(A) would imply that rhj = for all j > 1 since the distribution of m after any jump 
is unbiased. But this cannot be correct since then at long time differences, where the 
probability that at least one jump has taken place approaches one, the susceptibility x 
would have to drop to zero. Assumption (B) is also in general not correct; e.g. in the 
standard Bouchaud model the jump probabilities Pj clearly depend on the value of the 
starting magnetization since the distribution of the time until the first jump does. 

There is, however, one scenario in which these objections do not apply: if the exit 
rate r{E,m) from a trap of energy E and magnetization m is independent of h (and 
hence of m), both assumptions (A) and (B) are correct because the field biases neither 
the probability of not having jumped since the j-th jump, nor the probability of j jumps 
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having occurred. In the Bouchaud model, this field-independence of the exit rate holds 
when ( = K,/{1 + k). One then works out easily, by considering how the Gaussian 
distribution of m changes with every jump, that 

mj{...) = K^m^ + (3h{l - 0(1 + k){1 - K^) (100) 

The susceptibility is then x = /^(l ~ C)(l + ~ '^'') = ^ i^^)^ where the average 
is over the distribution of jumps / dE„ Pj{E^,t — t^)P{E^,t^). Since the correlation 
function is C = (/«■'), it follows that Tx = 1 — C and one has equilibrium FD behaviour. 
This is entirely consistent with the exact relation 

In the Barrat and Mezard model, assumption (A) can still be justified for k, ^ 1: 
one can show |12j that the exit rates are affected by the field via terms of 0{{1 — K,)h), 
which vanish in the limit of slow decorrelation. One is then allowed to find rhj 
from the distribution of m directly after the j-th jump, giving at low temperatures 
Thj{. . .) = n^m^ + 2h{l — K,^) 17 . However, for k, ^ 1 one needs a number of jumps 
j = (9(1/(1 — K,)) to see any significant response of the system, and similarly any 
significant decay of the correlation function. It is then not clear that assumption (B) 
can still be justified, since the small 0(1 — k) changes of each exit rate accumulated over 
this many jumps could still add up to a nontrivial m^-dependence of Pjit, h, E^, m„, t^). 

We note finally that, for T > 1/2 and k ^ 1 in the Barrat and Mezard model, one 
can state the argument for a straight-line FD plot without actually requiring assumption 
(B): histories with a finite number of jumps j do not contribute to x because for k —>■ 1 
the magnetization remains pinned to its initial value, whatever the field h. On the other 
hand, histories with diverging j always give the equilibrium susceptilibity because they 
pass through the equilibrated shallow traps at the top of the landscape. Because C is 
just the probability of j being finite as argued in Sec. 17.21 this gives the equilibrium FD 
relation x = Xeq(l — C) = j3{l — C). 

9. Conclusion 

We have studied trap models with observables that decorrelate slowly, by a factor k, ^ 1 
with each jump. Our motivation was to clarify and extend with an analytical study 
the interesting observations of Junier and Bertin for a spin model which leads to trap 
model dynamics of this type J7] . Our analysis showed that in the limit of long times, 
correlation and response functions are determined by two scaling functions: V{uj), which 
gives the scaling of the distribution across trap energies, and 7W(r), which depends on 
K and encapsulates the decay of the observable with scaled time. Both of these are 
determined by integro-differential equation which can be solved exactly at T = and 
by appropriate numerical iteration procedures otherwise. 

We focussed mostly on the Barrat and Mezard trap model which has entropic 
barriers that slow its dynamics at low temperature. In Sec. El we considered this first 
for fast decorrelation, k = 0, to extend previous results which had been limited to 
T = 0. The results already reveal changes in the dynamics at T = 1/2, i.e. half the 
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glass transition temperature Tg = 1 of the model: for T < 1/2 the aging behaviour is 
simple, with the average hopping rate scaling as the inverse of the age, while for T > 1/2 
a power law scaling obtains which is dominated by rare excursions of the system to the 
top of its energy landscape, i.e. to the shallow traps. Correspondingly the initial slope 
of the FD plot of susceptibility versus correlation is 1/Tes with Tcs = 1/2 for T < 1/2, 
and 1/T for larger temperatures. The long-time susceptibility x(^?^w) for t ^ i.e. 
for C{t,t^) — > 0, remains constant throughout the glassy region T < 1, effectively 
"freezing" to its value at T = 1 as temperature is lowered through the glass transition. 

In Sec. ini we looked at the complementary case of T = but slow decorrelation, 
K > 0. Here all scaling functions can be found exactly. The initial slope of the FD plot, 
1/Teff = 2, remains ^-independent, while the asymptotic (for C — > 0) slope 2k, grows 
linearly with k. In the limit k — 1 of slow decorrelation the two slopes coincide, and 
the FD plot becomes a straight line. 

Finally, in Sec. d we extended the analysis to general decorrelation (k > 0) and 
nonzero temperature. Scaling functions now need to be found numerically, but we were 
able to predict the exponents for their power-law decays in closed form. Surprisingly, 
these exponents depend on whether Glauber or Metropolis transition rates are used, 
contrary to the usual expectation that such choices only have minor quantitative effects. 
We were able to show explicitly that the initial slopes of the FD plots are k- independent, 
and that in the limit k — > 1 of slow decorrelation the FD plots become straight lines 
as argued qualitatively by Junier and Bertin [T7]. Their slopes, coinciding as they 
must with the K-independent initial slopes, correspond to FD relations with an effective 
temperature T^fi = 1/2 for T< 1/2 and equilibrium FDT for T > 1/2. We discussed 
in detail the approach to the limit /t — > 1; for T > 1/2 this is discontinuous because 
rare excursions to the top of the landscape achieve full decorrelation (and an associated 
equilibrium response) whenever n is below unity by a nonzero amount 1 — k. 

We generalized to the Bouchaud trap model, and more broadly any trap model 
with rates that are multiplicatively perturbed by the applied field, in Sec. |H1 Here it 
turns out that straight-line FD plots always have an equilibrium slope of 1/T. However, 
they do not necessarily arise even in the limit k —>■ 1, but rely on a specific assignment 
of the field-dependence of the transition rates which eliminates the effect of the field 
on the residence time in a given trap. This example illustrated that simple arguments 
for the existence of straight-line FD relations require some caution because the implicit 
assumptions can be difficult to justify. 

Looking at our results more broadly, it is intriguing to see that slowly decorrelating 
observables can produce non-trivial straight-line FD plots in trap model, something 
which can otherwise be achieved only by choosing rates that violate detailed balance jTU]. 
The significance of this is that such plots indicate that the effective temperature is (in 
the long-time, non-equilibrium regime) independent of the pair of observation times; 
this is one of the plausible requirements for Tggf to be a physically meaningful quantity. 
What remains unclear to us is why the notional "glass transition" temperature which 
slowly decorrelating observables seem to measure in trap models is half the usual 
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glass transition temperature Tg (= 1 in our units). One way of interpreting this 
result ^Tj is that slowly decorrelating observables only pick up non-equilibrium effects 
caused by entropic rather than energetic barriers. Alternatively, one could argue that 
slowly decorrelating observables detect only persistently slow dynamics. As soon as the 
system begins to exhibit intermittent behaviour, where episodes of fast dynamics in the 
effectively equilibrated parts of phase space alternate with periods of slow wandering 
around regions of low energy, effective equilibrium FD relations are recovered. Slowly 
decorrelating observables can thus help us to single out parts of the dynamics which 
can be meaningfully associated with an effective temperature, but by the same token 
can be blind to other aspects of the dynamics that remain clearly out of equilibrium. 
How this trade-off operates in other glassy systems is certainly worthy of further study. 
One challenge will be to integrate finite spatial dimensionality into the picture. One 
expects, for example, that the dynamics will then effectively decompose into that 
of independent subsystems whose finite size is set by a time-dependent correlation 
length [22] ■ Decorrelation can then not be arbitrarily slow: if each subsystem contains 
Nes degrees of freedom, say, one would expect k < 1 — 0{1/Ncs). The hmit of 
slow decorrelation we focussed on here would then require relatively large correlated 
subsystems. 
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